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Diseases emerge, persist and vanish in an ongoing battle for available hosts. Hosts, on the other 
hand, defend themselves by developing immunity that limits the ability of pathogens to reinfect 
them. We here explore a multi-disease system with emphasis on mutual exclusion. We demonstrate 
that such a system develops towards a steady state, where the spread of individual diseases self- 
organizes to a state close to that of critical percolation, without any global control mechanism or 
separation of time scale. For a broad range of introduction rates of new diseases, the likelihood of 
transmitting diseases remains approximately constant. 

PACS numbers: 89.75.Fb, 64.60.ah, 64.60.al, 05.65. +b 
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I. INTRODUCTION 

Many phenomena within materials science, physics, 
and biology are associated with percolation theory [TJ |2] . 
In particular, it has been shown that bond percolation is 
equivalent to the class of susceptible/infectious/recovered 
(SIR) epidemic models on a network [SHHj. 

In such an SIR model, all nodes start out susceptible 
to a new disease. If a node is infected, it will try to infect 
its neighbors for a fixed time r, after which it recovers 
and becomes immune to the disease. How widely each 
disease is spread on the network depends on the proba- 
bility p, with which a node infects each of its neighbors 
before it recovers. The probability p will, therefore, be 
directly given by the disease time r. If p is greater than a 
critical percolation threshold Pc, there is a finite probabil- 
ity that the disease will span the entire infinite network, 
thus becoming an epidemic. 

In nature, percolation phenomena are often found near 
the critical probability Pc [9 . A possible explanation 
of this is the concept of self-organized criticality (SOC), 
where complex systems drive themselves to critical states 
without the need for fine-tuning of the parameters [10]- 
[12]. Many models for self-organized percolation have 
been studied p!3HT7] . In these, the self-organization ei- 
ther arises as a result of very different time scales or 
through dynamics involving a global control mechanism. 
For instance, a percolation system can self-organize to the 
critical threshold by dynamically adjusting the probabil- 
ity p, such that the percolation cluster keeps growing at 
a specific rate [14 . However, this requires that all nodes 
on the network 'know' how fast the cluster is growing 
globally; a condition that is rarely fulfilled. 

In this paper we study an SIR model for the spread of 
multiple diseases that compete with each other. When 
many diseases are present, they may well influence each 
other by weakening of host immunity, or they may inhibit 
each other through cross-immunization or by mutual ex- 
clusion [ 181421] . Considering mutual exclusion only, we 
here show that the system self-organizes to a state close 
to the critical percolation threshold for a wide range of in- 
put parameters. That is, it exhibits self-organized quasi- 
criticality without any global control mechanism or sep- 
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Figure 1. At each time step, a random node 'coughs' on a 
random neighbor and thereby transmits a random one of its 
diseases. If the neighbor is not already immune to the disease, 
it is infected and becomes infectious for a time r. 



aration of time scale. 



II. MODEL 



In our model, diseases are spread on a 2-d square lattice 
with periodic boundary conditions and N = L'^ sites, 
each representing a host. At each time step the following 
actions take place: 

• With the small probability -^, a new disease origi- 
nates in a random node on the network. 

• A random node i and one of its four neighbors j 
are selected. If i carries any disease(s), a random 
disease is selected. If j is not already infected or 
immune to this disease, it is transmitted to j (see 

• After r sweeps over the lattice, j will be cured from 
the disease. 

The model uses a framework similar to the one recently 
developed in [22 , except that the present model allows 
each node to be infected by several diseases at the same 
time. Also, the present model has two input parameters; 
a is the introduction rate of new diseases on the network 



while the disease time r corresponds to the duration any 
node is infectious with a disease. 

A key element is that each node can only transmit one 
of its diseases at any given time step. Thus, a disease is 
less likely to spread from a node that carries many other 
diseases. The model can be run online as a Java- applet 
at cmol. nbi. dk/models/ disease/ MultipleDiseases. html. 

If a node is constantly infected with k diseases, the 
probability that it tries to infect a given neighbor with a 
given disease before it is cured, can be found to be 
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This probability corresponds to the percolation probabil- 
ity of a bond percolation system. When a ^ 0, diseases 
are rare and infected nodes will have only one disease. 
When r is large, this disease will have plenty of time to 
infect its neighbors, and will therefore spread in a circular 
manner with a broad rim of infected nodes and a solid in- 
terior of recovered nodes, very similar to the well-studied 
Eden growth [23-25 . When r is small, this disease will 
rarely manage to infect a neighbor before dying out. For 
r = 41n(2) ^ 2.77, we see from ([l]) that the probability of 
infecting a given neighbor is exactly Pc= \i which is the 
critical threshold of bond percolation on a 2d-lattice [26] . 
In this case, the spread of the disease becomes fractal- 
like and the regions of infected nodes are only one node 
thick. This behavior, which is shown in figure [2|a)-[2|c) 
for a system of size L = 256, is well-studied for single 
disease models [3]. 

When a > the analysis is made complicated by spa- 
tial and temporal variations in k and, therefore, p. For 
r > 2.77, the system will initially be supercritical, and 
all diseases will grow rapidly. Consequently, the average 
number of diseases per node (k) will increase and the av- 
erage probability (p) over the lattice will decrease. If (p) 
becomes less than ^, most new diseases will only spread 
to a handful of nodes before dying out. Thus, (k) will 
decrease and (p) will increase. This negative feedback 
mechanism will drive the system to a state, with (p) close 
to the critical percolation threshold, where disease sizes 
of all orders of magnitudes occur. In this state, the num- 
ber of diseases per node is close to Poisson distributed 
across the lattice, but with both spatial and temporal 
correlations in k. A high r will result in many diseases 
per node and a high a will make the system self-organize 
faster, but for a wide range of both parameters diseases 
will spread in fractal-like shapes, as can be seen in figure 

id)i;f). 

The average probability (p) of transmitting a disease 
can be measured directly by monitoring how many neigh- 
bors each node on average tries to infect with a disease, 
before it is cured. In figure |2|g)-|2|i) the development (p) 
is shown for various input parameters. It is seen that (p) 
converges to a value close to (p) = ^ . 




Figure 2. Spread of diseases on a lattice of size L = 256 
for different input parameters. Black nodes are healthy, and 
bright nodes carry many diseases. Accentuated clusters repre- 
sent a particular disease. In (a — c) only one disease is present 
{a ^ 0). U T ^ 2.77, diseases will spread in fractal shapes of 
critical percolation clusters. In the supercritical (r > 2.77) 
and subcritical (r < 2.77) cases, diseases will grow to span 
the entire network or quickly die out, respectively. In (d — /), 
a > the number of diseases per node will self-organize to a 
value where the transmitting probability is close to the crit- 
ical threshold of percolation p^ = ^. Thus, disease clusters 
have fractal shapes for a wide range of input parameters. The 
development of p is shown in (g — i). 



III. CRITICAL EXPONENTS 

To compare the model to a percolation system, the 
clusters of recovered and immune nodes were investigated 
for different sets of input parameters. For each disease, 
its cluster diameter, mass and exterior perimeter were 
measured. 

For critical percolation the cluster mass scales with the 
diameter giving a fractal dimension of D = ^ ^ 1.896 
[2 . In figure 3(a) it is seen that disease clusters have 
dimensions very close to this, for a wide range of input 
parameters. The fractal dimension is larger than || when 
- is small and vice versa. 

r 

The exterior perimeter is defined as the number of sites 
in the cluster that have one or more neighbors strictly 
outside the cluster. For critical percolation, it scales with 
the diameter with the critical exponent Dg = | [2 . In 
figure [3(b)] this is seen to be very close to the scaling of 
disease clusters. 

At the critical point, the cluster size distribution is ex- 
pected to fall off with the critical exponent — 1| ^ —1.05, 
such that small diseases occur more frequently than large. 
In figure |3(c) it is seen that the disease clusters fall off 
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(a) Cluster mass scales with the diameter with exponent close 
to the fractal dimension of critical percolation D ^ 1.896. 
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(b) External perimeter scales with the diameter with exponent 
close to that of critical percolation De = |. 
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(c) Cluster size distribution falls off with exponent broadly 
distributed around that of critical percolation —1.05. 

Figure 3. Critical exponents for disease clusters when a = 2 
and T = 5. Black lines show the exponents for critical per- 
colation, while the fitted values for the exponents are shown 
in the tables for different input parameters. The inserts show 
the experimental data normalized with respect to the critical 
percolation exponents. 




Figure 4. Steady state average transmitting probability (p) 
shown for different input parameters. When r > 2.77, the 
number of diseases per node self-organizes to a value, such 
that (p) is close to the critical probability of percolation. For 
low a, the system is slightly supercritical and, conversely, 
when a is high, the system is slightly subcritical. Fine-tuning 
of a parameter is necessary in order for the system to be truly 
critical. Note that the critical threshold is not necessarily at 
Pc = I due to correlations in the number of diseases per node. 



with exponents broadly distributed around this value, 
with steeper exponents when - is small. Here, the chance 
of getting a disease spanning the entire network is large, 
but the chance of a large disease suddenly dying out is 
low. 



IV. DISCUSSION 

In figure [4J the steady state probability (p) is shown 
as a function of the input parameters. The observations 
agree well with the characteristics of quasi-criticality [27] . 
The system self-organizes to a near-critical state, but a 
fine-tuning of a parameter (e.g. a) is necessary in or- 
der for the system to be truly critical. When a is too 
low, diseases are transmitted with a probability some- 
what larger that Pc^ and the system is supercritical - the 
disease clusters become 'heavy' with a fractal dimension 
above I^g = H^ an external perimeter dimension below 
D = ^ and high probability of forming an epidemic. 
When a is too high, the system is subcritical with low 
{p) and D, high Dg? and with low probability of forming 
an epidemic. 

It should be emphasized that, due to correlations in 
the number of diseases per node, the critical threshold 
of the multiple disease model is not necessarily equal to 
that of critical percolation Pc = ^ • Diseases will tend to 
'get stuck' and accumulate in regions where there are al- 
ready many diseases, while they will quickly 'sweep over' 
regions with few diseases. 

In the model, a node carrying k diseases will have a 



probability of ^ to pass on each of them. This couphng 
mechanism between diseases is not based on empirical ev- 
idence, but merely reflects that diseases compete against 
each other. However, the self-organization is robust to 
changes in the dynamics as long as diseases inhibit each 
other. For instance, if a node carrying k diseases has a 
probability of p- to pass on each of them, the case a = 1 
and r = 5 gives clusters with D = 1.875, De = 1.331, 
which indeed is close to the critical exponents of perco- 
lation. 



V. CONLUSION 

Using a recently developed framework for spread of 
many diseases [22 we have presented a simple multiple 
disease model that exhibits self-organized quasi-critical 
percolation [27]. The model is based only on local in- 
formation, having no global control mechanism or sepa- 
ration of time scale. Furthermore, the main feature of 
the self-organization is robust to changes in the extent 
to which multiple diseases inhibit each other. The basic 
mechanism employed in this model may be applicable to 
other systems, where "new" has an intrinsic advantage 
over "old", and where transmission capacity is limited. 
Thus the model may be equally valid as a rumor spread- 
ing model, where rumors compete for attention and be- 
come locally outdated. In this framework the predicted 
self-organization may help to explain the broad distribu- 
tions found in human social activities. 
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